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“What’s  the  buzz?”  -  Great  Promise: 

“Due  to  the  combination  of  a  high  speed  of  sound,  hard  optical  vibration 
modes  and  a  large  phonon  mean  free  path,  carbon  nanotubes  prove  to  be 
the  most  efficient  thermal  conductors.”  (D.  Tomanek,  2005) 


In  MD  simulations  of  Tomanek’s  group 
(2000),  the  peak  therm,  cond.  of  (10,10) 
SWNT  found  to  be  37000  W/m-K,  on  par 
with  the  highest  value  ever  observed  in  a 
solid  (diamond).  RT  value  of  6600 
W/m-K  is  also  very  high... 


“The  stiff  sp3  bonds,  resulting  in  a  high 
speed  of  sound,  make  monocrystalline 
diamond  one  of  the  best  thermal 
conductors.  An  unusually  high  thermal 
conductance  should  also  be  expected  in 
carbon  nanotubes, which  are  held  together 
by  even  stronger  sp2  bonds. . . 

The  rigidity  of  nanotubes,  combined  with  a 
virtual  absence  of  atomic  defects  or 
coupling  to  soft  phonon  modes  of  the 
embedding  medium,  should  make  these 
systems  very  good  candidates  for  efficient 
thermal  conductors”. 
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But  can  we  measure  and  utilize  it? 


•  “The  expected  high  thermal  conductivity  value  in  nanotubes  makes  its  direct 
experimental  observation  very  difficult.  Since  thermal  transport  is  likely  to  be 
dominated  by  phonon  scattering  in  the  contact  region,  the  most  important 
and  currently  unsolved  challenge  is  to  reproducibly  create  and  characterize 
thermal  contacts  to  a  nanotube.  Even  if  thermal  scattering  in  the  contact 
region  could  be  minimized,  the  net  thermal  transport  would  still  be  limited  by 
that  of  the  leads”.  (D.  Tomanek,  2005) 

•  “Interlayer  interactions  quench  the  thermal  conductivity  ...  by  nearly  1  order 
of  magnitude.  We  should  expect  a  similar  reduction  of  the  thermal 
conductivity  when  ...  nanotubes  form  a  bundle  or  rope,  become  nested  in 
multiwall  nanotubes,  or  interact  with  other  nanotubes  in  the  “nanotube  mat” 
of  “bucky  paper”.  (Berber,  Kwon  and  Tomanek,  2000) 


In  fact,  even  experimental  data  are  relatively  scarce 


Some  “standard”  references 


1.4 


(Hone  et  al,  1999).  Analysis  •§  o.6 
suggests  that  phonons  dominate  ^  04 
thermal  transport. 
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Experimental  Data  continued:  Morphology  and  Ordering  Matter 


23a 


Thermal  conductivity  of  the  “thick” 
annealed  sample  of  aligned  SWNT, 
measured  in  the  parallel  direction.  At 
300  K,  k  is  much  higher  than  in 
unaligned  material,  and  is  within  an 
order  of  magnitude  of  graphite  or 
diamond.  ( Hone  et  al,  2000) 
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Crossover  suggests  1-d  quantization  is 
indeed  up  to  higher  T  for  smaller 
diameters  ( Llaguno  et  al,  2002) 


Morphology  Effects  (Defects)  can  be  Improved  also  by  additional 
processing 


Anneanling  helps  by  eliminating 
impurities  and  healing  defects  ( Llaguno 
et  al,  2002) 
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Experimental  Data:  Multiwall  Nanotubes 


10  100 
Tempcralurc  (K) 

Individual  suspended  MWNT 
(d=14  nm)  per  Kim,  Shi, 
Majumdar,  and  McEuen  (2001). 
Maximum  k=3000  W/m-K  much 
larger  than  20  W/m-K  reported 
for  mats  (V7  et  al,  1999). 
Bundling  deteriorates  k. 
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k  ~  20  W/m-K  deduced  independent  of 
tube  length.  “Taking  the  volume-filling 
fraction  of  CNT’s  into  account,  the 
effective  thermal  conductivity  for  the 
MWNT’s  is  about  200  W/m  K”. 


Overall  sentiment  on  experimental  situation  seems  to  be 
that  room-temperature  values  of  k  vary  from  units  and 
tens  of  W/m-K  to  a  couple  of  hundred  of  W/m-K  for 
various  arrays,  the  results  ordinarily  improving  with  the 
decrease  of  the  bundle  size.  A  record  result  of  3000 
W/m-K  relates  to  individual  MWNT. 

Can  this  record  result  have  to  do  with  a  higher  density  of 
layers  of  the  “material”  of  the  sample  accompanied  by 
the  absence  of  the  need  of  entangled  3-d  heat  transfer? 

Also,  k  seems  to  be  an  increasing  function  of 
temperature  up  to  about  room  temperatures. 


•  Now,  what  about  our  understanding  of  phonons  in  NTs 
and  the  thermal  transport? 
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Model  vibrational  spectra  shown  using  spiral  symmetry 


Graphene  (10,10) 


Vibrational  spectra  of  individual 
SWNTs  are  basically  well 
understood: 


High-frequency  part  of  the  spectrum 
is  well  represented  by  zone-folding  of 
the  graphene  spectrum. 


Low-frequency  part  reflects  generic 
1-d  features  common  with  vibrations 
of  elastic  cylinders. 


Low-frequency  part  of  the  spectrum 
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features  4  gapless  acoustic  modes: 
longitudinal  and  twisting  with  linear 
dispersion  co  ~  k  and  degenerate  bending 
with  parabolic  dispersion  oo  ~  k2  as  well  as 
higher  energy  modes  with  quantized  gaps. 


In  the  “older”  picture  (however 
being  sometimes  used  even  now) 
bending  (transverse)  modes  were 
believed  to  have  oo  ~k. 
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Intertube  Phonon  Coupling  and  Dimensional  Crossovers 


Experimental  data  of  Hone  et  al  (2000)  on  specific 
heat  in  SWNT  ropes  suggests  that  intertube  phonon 


•Isolated  tube  phonons  cross  over  from  a  ID  regime 
where  only  acoustic  subbands  are  occupied  to  a  2D 
regime  as  higher  (optic)  subbands  are  populated. 

•  In  a  bundle  of  weakly  coupled  tubes,  phonons  are  3D  at 
very  low  temperature,  then  crossing  over  to  a  ID  regime. 
•If  the  coupling  were  strong,  the  ID  regime  would  be 
bypassed  and  a  quantized  1-d  phonon  spectrum  would  not 
be  observed. 


More  accurate  analysis  of  specific  heat  by  Popov  (2002)  does  not  seem  to 
contradict  the  conclusion  of  weak  intertube  phonon  coupling 


System 

Acoustic 

branch 

Phniinn 

dispersion 

Phonon  DOS 

Specific  heal 

Graphite 

LA.  TA 

Graphene 

LA.  TA 

D{ 

C^T* 

ZA 

Z>|»  =  cun  si 

or 

SWNT 

LA.  TW 

I >1  fkt )  =  enmt 

or 

TA 

Phonon  DOS  with  correct 
bending  modes 


6 


More  simulations  of  thermal  conductivity  of  SWNTs 


MD  simulations  of  Osman  and  Srivastava 
(2001)  also  (comp,  to  Tomanek’s)  reveal 
a  maximum  of  k  but  of  different 
magnitude  and  at  higher  T,  and  indicate  a 
dependence  on  tube  radius 


"Tern  para  lure-  IKh 


MD  simulations  of  Grujicic,  Cao  and  Roy  (2005) 


Results  are  closer  to  Tomanek’s 
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Despite  the  variations,  MD  results  show  a  similar  pattern  well  known 
in  the  theory  of  solids 


larger  crystal 


“Between  the  Umklapp  region  and  the 
boundary  scattering  region...  the 
conductivity  maximum,  whose  appearance 
is  perhaps  one  of  the  most  striking 
predictions  of  the  quantum  theory  of  solids”. 
(J.  Ziman,  1960) 

Qualitatively,  k  ~  C-v7 


At  high  T ,  it  is  phonon-phonon  interactions  (umklapp  scattering)  that  reduce 
mean-free  path  /. 

At  low  T,  umklap  processes  are  not  efficient  and  mean-free  path  is  overtaken 
by  the  sample  size  L.  T  -dependence  is  largely  that  of  heat  capacity  C. 


But  what  is  it  that  plays  the  role  of  L  in  MD  simulations,  where  sample  sizes 
are  at  best  hundreds  of  angstroms  with  or  without  periodic  boundaries? 
Estimates  for  mean-free  paths  are  normally  given  as  fractions  or  comparable 
to  a  micron.  Moreover,  one  can  make  simple  but  fundamental  estimates... 


Conductance  and  Conductivity 

Conductance  G  of  a  '’sample”:  Heat  current  J  =  -G  ■  AT 

iT  Conductivity  k  of  a  "material":  Heat  current  density  j  =  -K^— 

dx 

_  k  =  G^;  Length  L,  cross  -  section  S 

•  Additional  scatter  in  reported  results  is  due  to  various  choices  of  S 

Low-T  Ballistic  Conductance  Quantization 

(Rego  and  Kirczenow,  1998) 

In  the  absence  of  phonon  interactions  ( ballistic )  and  for  ideal  contacts,  each 
acoustic  mode  in  the  limit  of  low  T  brings  quantized  contribution  to  the  conductance 

G,  =^T  =  9.456 •  1(T13  (W/K2)  T;  G  =  N^G, 


3  h 


g  =  ^L  y 


In  general,  ballistic  conductance 


r  x  e  Tico 

dx - Z- - 

J  (px-W2  k„T 


a  -all  monotonic  segments  z 


kBT 
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Ballistic  conductivity  results:  A  likely  upper  limit 


Trying  to  rationalize  MD  results, 
consistency  check...: 

Knowing  vibrational  spectrum,  it  is 
very  easy  to  calculate  ballistic 
conductance,  no  other  scattering  but 
the  boundary  one.  Here  ball.  cond. 
is  shown  translated  into  k  for  a 
(10,10)  tube,  give  or  take  for  the 
accuracy  of  the  spectrum. 

Length  used  is  L= 1  micron  and 
cross-section  S= 250  angstrom2. 

To  get,  e.g.,  Tomanek’s  result  at 
r=100  K,  the  length  would  need  to 
be  increased  by  a  factor  of  ~  80.  So 
what  is  the  meaning  of  the  effective 
lengths  appearing  in  those  MD 
calculations? 


Temperature,  K 


Temperature,  K 


Interesting  comments  on  the  convergence  of  MD  results  upon 
the  increase  of  the  simulation  sample  size 
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"...  despite  the  fact  that  the 
phonon  mean  free  path  is 
considerably  larger  than  the  sizes 
of  the  computational  cells  used, 
an  apparent  convergence  in  the 
thermal  conductivity  can  be 
obtained”.  (Grujicic  et  al,  2005) 


“We  find  that  the  low-frequency 
vibrational  modes  of  the  lattice 
are  limited  by  the  size  of 
simulation  domain,  and  the 
thermal  conductance  of  an  infinite 
long  CNT  may  be  infinite”.  (Yao  et 
al,  2005) 


X  [W/ni  K] 


Intertube  Phonon  Coupling  and  Thermal  Conductivity 

Not  that  much  has  been  done. . . 


Berber  et  al  (2000)  base  their 
conjecture  on  the  analysis  of  the 
interlayer  effect  for  graphite: 

“Very  interesting  is  the  fact  that  once 
graphene  layers  are  stacked  in  graphite, 
the  interlayer  interactions  quench  the 
thermal  conductivity  of  this  system  by 
nearly  1  order  of  magnitude...  We  should 
expect  a  similar  reduction  of  the  thermal 
conductivity  when  a  nanotube  is  brought 
into  contact  with  other  systems”. 


Che,  Cagm  and  Goddard  (2000):  “We  also  carried  out  the  thermal  conductivity 
calculation  for  (10,  10)  nanotube  bundles  in  close  packing  condition.  The  simulations 
show  that  the  nanotube  bundle  has  very  high  thermal  conductivity  along  the  tube  axis, 
9.5  W/cm/K  [T=300  K],  which  is  comparable  to  simulated  graphite  in-plane  thermal 
conductivity,  10  W/cm/K.” 


Thermal  Conductivity  Enhancement  in  Composites 


MWNTs  in  Synthetic  Oil 
(Choi  et  al  2001) 
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“In-House”  (UTD  Nano-Tech)  Measurements 

{courtesy  of  Dr.  A.  Aliev) 


a  =0.2  W/m*K 


Polycarbonate  (Poly(Bispher  :il  A  carbonate)) 
[-(3C6H4A-C(CHJ),C#H4A-OCO-]n 
Ctystalline  thermoplastic,  tough,  stable,  transparent. 
Electrical  resistivity:  1017  Ohm*cm 
Dielectric  constant:  2.9-3.2. 

Thermal  cond.:  0.19-0.22  W/m*K(23°C. 

Density:  1.2  g/cm3. 


Concentration,  % 


Substantial  enhancement  in 
nearly  uniform  mixtures  of 
SWNTs  in  polycarbonate 


Much  more  modest  effect  for 
mixtures  in  conjugated  polymer 
MEH-PPV  (attributed  to 
agglomeration  of  SWNTs). 


Effective  Medium  Analysis  for  Composites 

{Nan  et  al,  2003,  2004 ) 


low  volume  fraction  /  result 

1  |  f'P  KdJSo. 

3  „+2at^ 

d  Km 

Ke  -  composite 
Km  -  matrix 
kc  -  nanotube 
p  -  aspect  ratio 
d  -  nanotube  diameter 
Kapitza  radius  aK  =  RKfcm 
due  to  interfacial  thermal 
resistance  RK 


i*-5n 
rf=15nin 


- —  d=25  HffL-  4 

- ^ 
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im 


-Some  exp.  data  can  be 
I  rationalized  (here  kc=3000 
and  Km=0.4  W/m-K) 


Extremely  important  role  of 
interface  resistance 


(b> 

d= 20  nm 

/f^SxtO^K/W 
— ji-  UXW 
~-p-  1500 

—  p=  2000  > 


0.0  0.2  0.4  0.0  0.0  1.0 

/m 


10  20  30  40  50 

JL  m2K/W> 


11 


Thermal  Interface  Resistance  between  SWNT  and  octane  liquid 

( Cahill  and  Keblinski’s  groups,  2003,  2004) 


Exponential  decay  of  a 
pulsed  temperature 
difference  leading  to  a 
close  estimate  for 
interface  resistance. 
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Steady  state  flow  between 
NT  and  a  sink  -  local  T  vs 
distance  from  NT. 

Dramatic  temperature 
gradient  at  the  interface  -  the 
largest  thermal  resistance  is 
with  the  nanotube-liquid 
interface  estimated  to  be  - 
3-1  O'8  m2K/W  for  NT  -  4  nm 
long  and  d=0.7  nm.  With 
liquid  k  ~  0.1  W/m-K,  the 
interface  resistance  is 
equivalent  to  having  NT 
surrounded  by  -  5  nm  thick 
extra  layer  of  octane! 


Vibrational  Mode  Analysis  of  Interface  Resistance 


WaveN  timber  k 

Spectral  temperature  of  nanotube 
bending  modes  as  a  function  of  mode 
number.  The  lowest  frequency  mode 
has  nearly  the  same  temperature  as  the 
surrounding  fluid  while  the  high 
frequency  modes  have  temperatures 
close  to  the  average  temperature  of  the 
nanotube. 


Decrease  of  the  interface  resistance  with  the  NT 
length  led  to  conjecture  that  it  is  coupling  of  low- 
frequency  bending  modes  of  NT  to  low- 
frequency  octane  vibr.  modes  that  provides  heat 
flow  -  all  interactions  being  driven  by  dispersion 
forces  (octane-octane  and  octane-NT). 

Spectral  temperature  analysis  confirms. 


Shenogin’s  et  al  (2004)  scenario: 

1.  The  heat  energy  first  flows  from  the 
high  frequency  modes  to  the  low 
frequency  transverse  vibration  modes 
(controlled  by  the  intrinsic  phonon- 
phonon  scattering  in  the  tube) 

2.  It  then  transfers  to  the  layer  of  octane 
liquid  adjacent  to  the  nanotube 
(controlled  by  the  coupling  between 
soft  modes  in  the  tube  and  in  the 
liquid) 
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Can  Interface  Resistance  be  Decreased  by  Increasing  Coupling  to 
the  Matrix?  -  Yes  but  only  by  so  much  (Shenogin  et  al  2004) 


Chemical  functionalization  -  introducing  chemical 
bonds  between  molecules  of  matrix  (octane)  and 
carbon  atoms  of  nanotube: 

1.  Interface  resistance  IS  decreased  7..  but 

2.  Defects  formed  on  NT  decrease  its  own  k 

Net  result:  A  limited  improvement  by  a  factor  of  2 
for  a  range  of  aspect  ratios  J 


Fraction  of  fruKlionilizcd  atoms 


Concluding  Remarks 

“Modest”  applications  of  nanotubes  for  thermal  management  (like  100- 
200%  increase  for  composites  with  low  conductivity  matrices)  look  quite 
feasible. 

However,  in  order  to  evaluate  prospects  of  more  aggressive  goals  fully 
exploiting  the  claimed  potential,  much  more  research  work  is  needed  to 
provide  a  firm  understanding  of  issues  involved.  Among  the  problems 
could  be: 

•  Details  of  the  physics  of  individual  tube  behavior,  especially  the  low-T 
regime 

•  Microscopic  picture  of  intertube/interlayer  interactions  and  scattering 

•  Microscopic  study  of  contacts  with  various  media  (leads) 

•  Possibility  of  (self-)  assembly  of  contacts  and  environments  that  would 
be  beneficial  to  the  therm,  conductivity 

•  Optimization  issues  given  the  understanding  achieved 
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